Generating Function for Particle-Number Probability Distribution in Directed 

Percolation 



Lucian AntorQ 

Department of Physics and Astronomy, University of Manchester, MIS 9PL, U.K. and 
Institute of Atomic Physics, INFLPR, Lab 22, PO Box MG-36 R 76900, Bucharest, Romania 

Hyunggyu Parl^l and Su-Chan ParlJl 

School of Physics, Korea Institute for Advanced Study, Seoul 130-722, Korea 
(Dated: February 2, 2008) 

We derive a generic expression for the generating function (GF) of the particle-number probability 
distribution (PNPD) for a simple reaction diffusion model that belongs to the directed percolation 
universality class. Starting with a single particle on a lattice, we show that the GF of the PNPD can 
be written as an infinite series of cumulants taken at zero momentum. This series can be summed 
up into a complete form at the level of a mean-field approximation. Using the renormalization 
group techniques, we determine logarithmic corrections for the GF at the upper critical dimension. 
We also find the critical scaling form for the PNPD and check its universality numerically in one 
dimension. The critical scaling function is found to be universal up to two non-universal metric 
factors. 
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I. INTRODUCTION 

In the last decades, the statistical mechanics of non- 
equilibrium models characterized by large scale space- 
time fluctuations and absorbing phase transition has 
been studied extensively with numerical and analytical 
tools . The directed percolation (DP) universality class 
covers a large variety of models of this type, which have 
short range interactions, a single absorbing state, and a 
scalar order parameter 0, Q- The DP conjecture has 
been extended to models with multiple components and 
even multiple absorbing states, unless there is an addi- 
tional symmetry and a conservation law. Nevertheless 
we should mention that there is at least one exception in 
the case of more elaborated interactions 0, IE ■ 

Besides various numerical techniques, the response 
functional formalism combined with the renormalization 
group (RG) techniques is one of the powerful analytical 
tool to explore the asymptotic properties near criticality 
for these models 0H H IH E3- The DP universality 
is described by the Reggeon field theory which has been 
investigated in details by many researchers pj. Its up- 
per critical dimension is known to be 4 and the values 
of the critical exponents were obtained up to 0(e 2 ) with 
space dimension d = 4 — e 0, 0] • The leading logarith- 
mic corrections at the up per critical dimension are also 
available in the literature jl3l . ITl ] and tested numerically 
[l5j . Recently the same field-theoretical formalism was 
successfully applied to identify and calculate the survival 
probability and the mean number of particles, which are 
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the central quantities of interest in the damage-spreading 
(DS) type dynamics E^. 

The aim of this paper is to explore in the same frame- 
work a more general and useful quantity in the DS type 
dynamics; the generating function (GF) of the particle- 
number probability distribution (PNPD). In principle, 
all moments of the number of particles including the sur- 
vival probability can be derived from the GF. We start 
from a lattice description of a reaction diffusion model 
that belongs to the DP universality class. Then, we use 
the bosonic operator formalism to represent the mas- 
ter equation and the associated field theory to define the 
generating function (GF) of the PNPD. With a single 
particle initially, we derive the GF of the PNPD as an 
infinite series of multi-point Green functions. This series 
can be summed up explicitly under the tree (mean-field) 
approximations into a complete form. 

Assuming the generic critical scaling property of the 
Green functions, we obtain the critical scaling form of 
the PNPD. Using the RG techniques, we calculate the 
leading logarithmic corrections for the critical GF at the 
upper critical dimension. Our results are consistent with 
recent reports for the survival probability and the average 
number of particles by Janssen and Stenull [l4l fill . 

The conventional RG analysis suggests that the critical 
scaling function for the PNPD should be universal up to 
two non-universal metric factors which microscopic de- 
tails of the models are integrated into. We check this 
universality numerically in one dimension for two vari- 
ants of bosonic DP models and also for the hard core DP 
models where the site occupancy of particles is restricted 
to one. Numerical simulations show that all the models 
form a universality very nicely at criticality with only two 
adjusting parameters. 

The paper is organized as follows: Section [H] intro- 
duces the model and defines the GF for the PNPD. Sc- 
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cion IIIII is devoted to the derivation of the series expan- 
sion of the GF in terms of connected Green functions. 
In Sec. IIVI we derive the mean field solution and present 
the logarithmic corrections to the scaling in the critical 
state at the upper critical dimension for the first three 
moments of the PNPD. In Sec. we analyze the scal- 
ing form of the critical PNPD, using the RG equations 
and test it numerically in one dimension for both bosonic 
models and hard core models. We present the conclusion 
in Sec. I VII and the technical details for the calculation of 
the logarithmic corrections in the Appendix. 



for a given initial state |^(0)) 0,@jI3j where |0) repre- 
sents the vacuum state and O is the operator represen- 
tation of the quantity O. Similarly, one can easily derive 
the probability that the system has n particles on the 
lattice at time t as 

1 / \ n 

P(n,t) = ^<0|(I>) e-* ff |*(0)) , (5) 

i 

and its generating function (GF) is denned as 

^( S ,i) = ^ S >(n,t)-(0|e£. sa 'e- tH |vI/(0)) . (6) 

n=0 



II. GENERAL FORMALISM 

We consider a reaction diffusion model (the Gribov 
process) 0, that belongs to the DP universality. It is 
defined on a d-dimcnsional lattice with N d sites which 
may be either vacant or occupied by particles. The evo- 
lution rule is summarized in the stoichiometric notation 
as 
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> 0A , 




A 


A a 






2A 


Ac 


A, 


(1) 


A 


A b 


2A, 





where Di is the nearest neighbor hopping (diffusion) rate, 
A a is the annihilation rate for a particle, A c is the co- 
agulation rate for two particles at a site, and Ab is the 
branching rate. 

A microscopic state of the system is defined by the set 
of particle numbers at each lattice site: n = {n^}. The 
time evolution of the probability distribution P(n, t) can 
be described by the master equation for the state vector 

!*(*)) = En F ( n > *)l n ) at time t: 



4g) 
dt 



-mm 



(2) 



The evolution operator H can be written in terms of 
bosonic operators as 

H = Di y^(a| - a\){ai - Qj) - A a y^(a, - a]a t ) 

-A c ^(o|a, 2 - a^a?) - A fc ^(a^a 4 - a+iOiX?) 



where ai(a|) is the annihilation (creation) operator of 
bosonic particles, satisfying the commutation relation 
[otjo]] = Sij. 

The average of a quantity O at time t is computed by 
the formula 



(O)t = <0|eS*°*&r tH |*(0)) 



III. THE EXPANSION OF THE GENERATING 
FUNCTION 

It is well known that the calculation of the average 
becomes simpler by moving the factor { ai in Eq. Q 
to the rightmost position. Since & a (V = (a) + l)e a , this 
move results in shifting — * a' + 1 in the operators O 
and H . We use the same trick to calculate the GF in 
Eq. JBJ), and then 

F(s, t) = (0\e^^ s ^ a - e- tH ^e^^^ |*(0)) , (7) 

where the shifted Hamiltonian H s (a', a) = H{a* + 1, a). 

Through the standard coherent-state path integral for- 
mulation one can find the shifted action S(cj),(j)) 
such that the path integral is weighed by e~ s as 



S@,<t>)= I dt J2& d t<t>i + H s $,< 

(a) 

{( A ° - Xb)fa<f>i+Xcfa<f>i-h$<l>i+K$4>i} 



dt 



(8) 

For the calculation of the asymptotic properties, the ac- 
tion is taken formally in the continuum limit with the 
transformations 



E- 

i 

<f>i{t) 



j fxi~\-be 

(x,t), ^{t) 



Hx,t) , (9) 



D = b 2 D u fi = \ a - A h , g' = b d \ c , g" = X b , 

where b is the lattice constant and ej is the unit vector 
along the j axis. In the continuum limit, the action reads 



S 



dxdt [(j)(x, t)dttp(x, t) + D(Vcj)(x, t))(V</>(a;, t)) 



+(j,<t>(x, t)(j>{x, t) + -g'<f>(x, t)4> 2 (x, t) 



g"cj) 2 (x,t)(j)(x,t) 



(4) 



(10) 



3 



where the term g'^tf) 2 is dropped as being irrelevant. 
The coupling constants g' and g" can be reduced to one 
coupling constant g = y 'g' ' g" by rescaling the fie lds as 
follows: 4> ~ * a 4>i 4> ~ * a ~ 1 <fii with a = y/ g" / g' ■ The 
rescaled "symmetric" action is then invariant under the 
time reversal transformation: t — ► —t, (f> — ► —<f>, 4> — > — <^>. 

Consider the GF with one particle initially: 1^(0)) = 
aj|0). Then, we have 



F(s,t) 



,£i(«-i)«<p-«.f„t 



1 + (0| 



,E<(a-l)«i p -*».„t 



(* - 1) P 



E 



(oS + l)|0) 

Slo> 



5(o)|o> 
(ii) 



where we used the relations H s \0) = and (0|if s =0 be- 
cause each term in H s has to the left and a to the 
right. In the path- integral formulation, the above cor- 
relation functions are simply connected Green functions 
(cumulants) with the shifted action in Eq. (JHJ). In the 
continuum limit with the symmetric action, the GF is 
written as 



F(s 



*) = ! + £ 



(g - 1) P P 



OL 



x dx 1 ,...dx p G {p ' 1 \x 1 ,t,...x pi t;0,0) , (12) 



where 

G (p ' 1} (z x , i, . . . x p> t; 0, 0) = (<f>(xi ,*)... 0(s p , i)0(0, 0)) 

(13) 

are the connected Green functions for the symmetric ac- 
tion and the factor a = \J g" I g' stems from the rescaling 
of the fields. 

Close to the critical point, we can use the renormalized 
Green functions of which the scaling behavior is known. 
Then, we have 

G^ p ' 1 \xi,t, ...x p , t; 0, 0) 

= zfzfG<i' 1 \x u t,---xp,t;0,o) 



p/2 7 l/2 + -P^-^ n (p,l) 



zy'zy t 



n "\\g^{xx/t 



l/z 



■ Xp/t 1 /*) , (14) 



where Zq, Zi are the fields renormalization factors that 
are equal for the time reversal symmetric theory and we 
have used the scaling property of the renormalized Green 
functions close to the critical point as 

G^ 1 \x 1 ,t 1 ,---Xp,tp ] x 0l t Q ) = W P+P ' 
x Hdal^x^lfiph,...) , (15) 

where (i is the deviation from the critical point, and /3, /?', 
vu, u±_ and z are the standard critical exponents defined 
in Ref. Q. 



a) 




c) 



b) 



FIG. 1: The first three tree diagrams used in the mean-field 
approximation. 



The generating function can be written now as 
F(s,t) = F(s,t) - 1 



( S - 1)P aP- 1 Z p J 2 zl'*t- p{ 



E 

p=i *" 

C^zT-z 1 ' 2 



■ A (P,1) 



j2 {s ~J )P (azi/y-H Ai -^AM 

(16) 



.p=i 



where 



A^ = dx 1 ...,dx p g(P< 1 \x 1 ,...,x p ) . (17) 



This expression shows that F(s,t) can be renormalized 

1/2 1/2 

by multiplication with Z^ Zj after the parameter 
a is renormalized as ccr = zV 2 a . 



IV. MEAN FIELD APPROXIMATION AND 
LOGARITHMIC CORRECTIONS 

A. Mean-field solution 

In Eq. (|16f) . the GF is expressed as a formal series 
whose convergence property depends on the coefficients 
A^P' 1 ', As a first step to study this series, we calculate 
the mean field (MF) solution for the GF. In this approxi- 
mation, we consider only the tree diagrams for the multi- 
point correlation function G^' 1 ) in Eq. I|12fl . 

The lowest order approximation of the (p, 1) Green 
function can be obtained by connecting p — 1 vertices of 
type (j) 2 4> with free two point Green functions, as shown 
in Fig. n where each line is associated with the free two- 
point (1, 1) Green function in the momentum space as 



G(k,t) = 6{t)e- t(Dk2+ti K 



(18) 



For the sum of tree diagrams we can easily obtain an 
integral recurrence equation for the GF by resuming the 
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diagrams after the first vertex from the right in the tree, 
which is 

F(s, t) = F{s, 0)e~ flt +-g" / dt'F(s, t - t') 2 G{k = 0, t'), 
2 Jo 

(19) 

that is equivalent with the differential equation 

±F( S ,t) = -vF( S ,t) + ±g"F( S ,tf. (20) 

With the single particle initial condition F(s, 0) = s 
(F(s, 0) = s — 1), one can find the exact solution for the 
GF as 



F{s,t) = 1 + 



(s - l)e-^ 



l-( s -l)|-(l- e -Mt) 



(21) 



By inverting the GF in Eq. ©, we find the PNPD as 



p(n,t) 



-fti 



(l + |_(l_ e -^))2 yi + |_(i_ e -Mt) 



2/j ( 



(22) 



for ft > 1. At criticality (/x = 0), the PNPD behaves in 
the asymptotic limit (large i) as 



(23) 



The survival probability p s (i) that the system is active 
at time t can be derived from the expression p s (t) = 
1 - p(0, i)= 1 - F(0, t). From Eq. JHJ), we find 



i>.(t) 



l+|_(l_ e -Mt) 



(24) 



We also find the mean number of particles N(t) = (n) = 
Y,n n P( n ^) = dF/ds\ s= i as 



N{t) = er^. 



(25) 



Similarly, all higher moments (n k ) can be obtained by 
differentiating the GF with respect to s and taking the 
s — > 1 limit. 

At criticality, the survival probability decays alge- 
braically as 



1 2 -i 

Ps[t) ~ l + g"t/2 W J' 1 ' 



(26) 



while N(t) remains a constant of unity. This is consis- 
tent with other MF predictions for the dynamic scaling 
exponents: 5mf = 1 and T)mf = 0. 

Near criticality, we find, for the large t limit, that 
p s (t) ~ (2fi/g")e^ ,lt in the absorbing side (p, > 0) and 
p s (t) w (2\p\/g")(l + e-M*) in the active side (p, < 0). 
These results are fully consistent with the recent results 
by Janssen [l6[ . 



B. Logarithmic corrections to critical scaling 

At the critical dimension, where the coupling constant 
is dimensionlcss, the MF solution l|21|) can be improved 
systematically at large time with the help of the renor- 
malization group equations by adding the one-loop cor- 
rection. A detailed account of this procedure has been 
presented recently in reference |l4l| . In our case we use 
one-loop approximation for the GF, combined with the 
two-loop Wilson functions taken from literature. The 
details of the calculation for one-loop GF and its asymp- 
totic scaling derived with the help of RG equations are 
presented in the Appendix. In this section we present, for 
the critical state, in a parametric form, the leading order 
expressions of the survival probability, p s (t), and of the 
first two moments of the particle number distribution, 
N(t), -/V 2 (t), which might be of interest for numerical ex- 
periments. 

From Eq. [jA~23|) we have at large t 



Ps{t) 



N(t) 



N 2 (t) 



-F(0,t) 
4N' 

0(w- 5 / 6 t- 2 ) , 
dF(s,t) 



W - I e cow2+O{wi) (1 + (31nA - 9/2)u; 2 ) 

(27) 



s=l 



N'w~ 1/3 e Clw2+0( - wl \l + (lnX + l)w 2 ) , (28) 
d 2 F(s, t) 



N(t) 



1 



x (1 - (lnX - A)w 2 ) + Oiw- 1 / 3 ) 



(29) 



where w depends on time t through Eqs. (|A.25I IA.26fl . 

In the above expressions, the universal constants Co, ci, 
C2 arc defined in terms of the coefficients of the Wilson 
functions, Eqs (|A.7t IA.9|l . and have the values 



c = l(7-341ni 



1 / 4 
cx = -— [ 25 + 322 In - 



144 

1 
144 



1 / 4 
c 2 = -— [ 71 + 542 In - 



(30) 
(31) 
(32) 



The remaining N' , D, X$ and to are non- universal con- 
stants whose detailed origins can be explored in the Ap- 
pendix. 

We make the observation that Eqs. <|27I28|I are equiv- 
alent with the expressions of Ref. after the change 



of variable w 



(w/8) 



1/2 
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V. CRITICAL SCALING FUNCTIONS FOR 
THE PNPD 

Since the GF of the PNPD is renormalizable, see 
Eq. (|16[) . we can use the RG equation to find its scal- 
ing properties in the critical state. At large t, we can 
extract the scaling property of the PNPD through the 
GF scaling behavior. 

A dimensional analysis with Eq. (|16|) suggests the fol- 
lowing form for the GF 



F(s,t) =yf(D K \ya K - d / 2 ,u) 



(33) 



where is an arbitrary length, u is the dimensionless 
coupling constant (see Appendix), and y = s — 1. 

The independence of the bare F on k yields the renor- 
malization group equation 



K lL (y z <i>f( DK2t >y aK ~ 

or, in more details, we have 



d/2 



,u)\ =0 



(34) 



x yf(DK 2 t,yaK- d/2 ,u) = 



(35) 



where 70 = nd K \nZ^ } j D = K.d K \nZ D , -f a = 7^/2 and 
fi(u) = nd K u, with Z^, Zd the standard renormalization 
factors Q. 

Equation l|33|l gives the following relations for the dif- 
ferential operators 



d d d d 

ok at 2 da 



D 



a 



3D 
d d 

a ^r = y~a — 1 ' 

aa ay 

hence the generic solution at the fixed point is 

1% ( d / 2 — '$/ 2 \ 



(36) 



(37) 



where starred parameters are values of corresponding 7 
functions at the fixed point. 

Using the standard RG technique the function / 
up to 0(e) and a multiplicative non-universal amplitude 
can be obtained by replacing in Eq. I|A.6(I u 2 — > e/12, rt 
-► -C* y yi (d/2 - 7 * /2)/(2+7 ° ) and t/t -> 1; where C v > is 
a second non-universal amplitude, 7^ = — e/6 and 7^ = 
-e/12 up to 0(e) 

Since the difference between p(n, t) and p(n + 1, t) for 
large t is small, the generating function F can be ap- 
proximated by a Laplace transform with the parameter 
lo = —Ins. Applying the inverse Laplace transform in 
which we approximate y f» — lj, we obtain 



p(n,t) = t- d ^g 



t d/(2z)+, 1 /2 



(38) 



o 




FIG. 2: Critical scaling function for the PNPD in one di- 
mension with full lines for boson models (pair annihilation 
and pair coagulation) and with dashed lines for hard-core 
particles models (the contact process and the bond directed 
percolation). The curve corresponding to the contact process 
is chosen for reference, and the other data sets are plotted 
for the best overlap with the reference curve by varying the 
constants Ci and C2. 



where we have introduced the exponents z = 2 + 7^, and 
r/ = —7^/ (2+7|j) defined by scaling properties of the two 
point Green function Equation 138|) along with the 
relation p s (t) = J dnp(n,t) gives the survival exponent 
8 = (d/z — T])/2, which is the well-known hyperscaling 
relation. 

The scaling analysis predicts that the critical PNPD 
at large t has the form 



p{n,t) 



c 2 



g[Ci 



t ri+8 



(39) 



where C\ , C2 are two non- universal metric factors which 
depend on the values of the microscopic parameters and 
g is the universal scaling function. The two-scale factor 
universality is worth while to be mentioned at this point. 
Since we are interested in the dynamic behavior of the 
system, there are, in general, three independent metric 
factors which are related to the number of particles, the 
time, and the distance from the critical point, respec- 
tively. The third metric factor does not appear in Eq. 
(|2HJ), because the system is at the critical point. This 
should be com par ed to the two-scale factor universality 
in equilibrium [l9j . 

In one dimension, it is relatively easy to test numer- 
ically the above scaling form. We performed Monte 
Carlo simulations for the bosonic model |2(| described 
by Eq. Jy| at two points in the parameter space 
(Di, A a , Afc, A c ) on the critical manifold: Di = 0.1 and 
0.2, A a = (1 - 2A)p 0J h = (1 - 2£>j)(l -Pc), A c = 0.5 
with p c = 0.11056(1) and 0.13345(4), respectively. We 
measured the PNPD at two different times t = 4000 
and 8000 for 2 x 10 7 samples. As shown in Fig. [21 wc 
overlap the scaling functions almost perfectly by varying 



the constants C\ and C 2 , which implies that the scaling 
function g is universal. The numerical values of the ex- 
ponents 77 and 8 are taken from Ref. [2l]]: ">] = 0.313686 
and 8 = 0.159464. In the same figure, we also plotted 
the scaling function of a variant model which has a pair 
annihilation process 2A — > instead of the pair coagu- 
lation 2A — > A. This one also overlaps perfectly with the 
others. 

As pointed out in Introduction, there are many other 
models belonging to the DP universality class. In the 
field theoretical formalism, they differ by extra irrele- 
vant coupling terms in the action besides the terms in 
Eq. i|10|) . Among these models, of a particular interest 
are the hard-core models where the dynamic rules accept 
only one or zero particle per site. As such examples, we 
have collected the PNPD data for the contact process 
and the directed bond percolation 0| in one dimension 
at the same times as for the previous data. Their scaling 
functions can be also overlapped very well with those of 
the bosonic models in the same manner. We have noticed 
a minute deviation between the two scaling functions at 
small values of the scaling variable Cin/t v+s . but a more 
detailed analysis shows that this deviation is due to cor- 
rections to scaling which decay very slowly in this region. 



VI. CONCLUSIONS 

We have shown that the generating function of the 
particle-number probability distribution of a DP model 
can be expressed as a series of cumulants of the associated 
field theory. Using the recurrence properties of the tree 
diagrams, we have calculated explicitly the mean field 
solution and the logarithmic corrections for the PNPD 
at the upper critical dimension. 

From the renormalization group analysis, we have 
shown that the scaling function for the PNPD is uni- 
versal in the critical state and the microscopic details of 
the model are included in two non-universal metric fac- 
tors. We have tested the universality numerically in one 
dimension for both boson and hard-core particle mod- 
els. Numerical simulations show that all models form 
the universality with only two adjusting parameters in 
the critical state. 
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FIG. 3: A generic one-loop diagram. The thick line represents 
the MF solution of the GF whose final time is t. 



APPENDIX: ONE-LOOP CALCULATION AND 
LOGARITHMIC CORRECTIONS 

The logarithmic corrections to the critical scaling at 
the upper critical dimension can be obtained using the 
RG equation (|35|l . For the calculation of logarithmic cor- 
rection of the critical GF up to the sub-leading order, we 
need the expression of the critical GF up to one-loop or- 
der and the Wilson functions up to two-loop order. The 
Wilson functions have been calculated previously in Ref. 
[Tlj . hence we have only to calculate one- loop correction 
to the critical GF. 

The generic topological structure of a one-loop dia- 
gram is shown in Fig. [3] To sum up all such diagrams, 
we need the onc-particle Green function with insertions 
of trees whose end points are located at the same final 
time t. The graphical representation of such calculations 
is shown in Fig. 

At zero mass, the one-particle Green function with in- 
sertions of p trees is 



K p (t 2 ,ti,k]t) = 
{ag) p / dri \ dr 2 

Jti Jti 

= -Dfc 2 (t 2 -t!) 



'2 



d,T p T(t - n) . . . T(t - t p ) 



x e 



2'' 



log 



l + r(t-t 2 ) 



-r>fc 2 (t 2 -*i) 



l + r(i-ii), 
where t 2 > t\ is assumed and 

T(t) = (s-l)/(l + rt) 



(A.l) 
(A.2) 



is the MF solution of the critical GF, i.e., F(s, t) at fi = 
and r = — (s — l)ag/2. Hence the sum over all tree 
insertions is 

K(t 2 ,t-L,k;t) = ^K p (t 2 ,tx,k;t) 

p=0 



l + r{t-t 2 ) \ 2 e _ Dk 2 (t2 _ ti) 
1 + r(t - ti) 



(A.3) 
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With the above Green function the sum over all one- 
loop diagrams has the form (see also Fig. [SJ. 

x K 2 (t 2 ,t 1 ,0;t)K(t ll 0,0;t)(t 2 -t 1 )- d ^ 2 , (A.4) 

where the integral over momentum was done. The in- 
tegrals in the foregoing expression are standard and we 
have, using dimensional rcgularization, 




L = — 



u 2 {s-l) 
(1+rt) 2 

2 rt ft 
'el + e/2 \i 
1 



1 



2 ft 



l-e/2 \e \i 
* _ 1 



rt 



2 + e/2 



l + e/2 
2rt - 21og(l + rt) - 



(rt - log(l + rt)) 
r 2 i 2 



1 + rt 



O(e) 



(A.5) 



where e = 4 - u = \fK±(g/ D)K~ t/2 /A with A" 4 = 
l/(87r 2 ) is the dimensionless coupling constant, and 
t^ 1 = 8ttDk 2 with k an arbitrary length scale. The 
renormalized GF can be obtained using the standard 
renormalization constants = 1 + (2/e)u 2 and Z g = 
1 + (5/e)u 2 ; see Eq. Ijlfijl . Hence, the renormalized GF 
up to one-loop order is 



F(s,t) = 



1 + rt 



1 + rt 



1 + log(t/t) 



-2 log(l + rt) + rt ( 3 log(t/*) - ^t^T 

1 + ?'t 



(A.6) 



FIG. 5: The diagram of one-loop corrections to the GF. 
The thick line whose terminal point is t is the MF GF and 
the double-line with arrow is the one-particle Green function 
with insertions of all trees. 



with the observation that in Ref. the Wilson func- 
tions are written for the square of our coupling constant 
multiplied by a different factor. 



The characteristic system to be solved is (see Eqs. (j35 



dp 
da(p) 
dp 

dt{p) 
dp 



(-2 + ^(u)/2)a(p), 
(2 + lD (u))t(p) . 



(A.10) 
(A.11) 
(A.12) 



For the calculation of the logarithmic correction we use 
the Wilson functions in two-loop order from Ref. 01 



0(u) = 6u 3 - - (l 69 + 106 In ^ ) u 5 



1d(u) 



-2u 



2 , 2 (6-91n-^- 4 



i(i7- ata |). 











"t 


il 




- 




+ n 






tl 







(A.7) 
(A.8) 
(A.9) 



K (4 2 , h,k;t) = G(k,t 2 -ti) + Ki(t 2 ,ti,k;t) + K 2 (h,t 1 ,k;t) +■■■ 

FIG. 4: The one-particle Green function with insertions of 
trees. 



Using dp/ p = dw/P(w) 1 we have 



t(w) = tC t w a " cxp[-l/(/? 3W 2 ) + b Q w 2 + 0(w 4 )} , 

(A.13) 

a(w) = C a w ai exp[l/(/3 3 w 2 ) + h x w 2 + 0(w 4 )} , (A.14) 
^Q-dw' = N'w a2 exp[fo 2 w 2 + 0(w 4 )} , (A.15) 



where Ct, C a are determined from the initial condition 
t(w = u) =t, a(w = u) = a. We denote by 4>i and Di 
the coefficients in the expansion of the Wilson functions 
P( w )-> lc/>(w) and 7_d(w) that generically can be written 
Q = J2i c QA ul > similar to the notations used in Ref. [l4j . 
The constants ao, Oi, a 2 , bo, b\. b 2 arc defined by the 



8 



following relations 
1 



ao = -^(^ 2 /33 - 2/3 5 ) , 
a i = i(02/3 3 + 4/3 5 ) , 

^ (^4/3 3 - D 2 (3 5 ) 



bo 
bi 



1 



2^2(04/33 - 02/3 5 ) ■ 



pi ' 

3 



(A.16) 

(A.17) 
(A.18) 
(A.19) 

(A.20) 
(A.21) 



Under scaling the generating function transforms as fol- 
lows: 



F(s, t; w, a) = exp 



dw'^(w')/P(w'} 



x F(s, t(w); w, a(w)) , 



(A.22) 



where we can eliminate w with the condition t(w) = X^t 
with Xq an arbitrary constant of order one. The GF at 
large t is 



F(s,t) 



N'(s - \) w ^e b2w2+0 ^\l + QnX + l)w 2 



S(w, s, t) 

-(1-S(w,s,t)) ( 31ogA - 



S(w, s, t) 



l + logA + 21og(S'(w,s,t)) 



!) 



2 2S(w,s,t) 



where 



S(w,s,t) = 1 - -(s- l)Dtw ao+ai+1 e^ bo+b ^ w2+ °^ , 



D = CaCfDK^ X ^ 2 and for large t we have from Eq. (|A.l 



—i 11 
w = t ' exp 



^i^I + 0(r- 2 ln 2 r) 
4 t y ' 



/3 3 Ht/t ) 



(A.23) 
(A.24) 

(A.25) 
(A.26) 



with t = {X /C t )t. 
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